Auto-Forecast
This module performs a DCA autofit on all streams (oil, gas and water) of selected wells in a given project. You can have multiple versions of these DCA autoforecasts (for different wells or different ranges of autofit parameters) within a project.
How does the autofit algorithm work?
The autofit algorithm finds the timestep with the maximum production data and performs a single-segment or multi-segment best fit from that timestep and rest of historical time.
Here's how to create a new autoforecast case:
- Click Auto-Forecast in the navigation panel.
- To the upper right, click ADD AUTO-FORECAST.
- Provide a name and select all the wells by clicking the check box at the top of the well list (left of Well Name).
- Click SAVE to the lower right.
You can store multiple auto-forecasts per project
This main autoforecast page gives you an overview of the created auto-forecasts, when they were created, who created them and what wells went into creating them.
After the autoforecast is created, there are three main steps: 1. Initialize, 2. Autofit and 3. Summary.
1. Initialize
This is where the autofit is initialized with both fit ranges and other relevant assumptions.
1.1. b-factor switch
The switch at the top of the screen dictates what should happen if the b-factor min or max is the converged solution during the autofit.
- When switch is ON: if the best fit solution hits the upper, or lower, b-factor bound, do a best fit with the default b-factor is done instead.
- When switch is OFF: if best fit solution hits the upper, or lower, b-factor bound, that is what is used.
1.2 Autofit Decline Type
You can use the Autofit Type dropdown to select between the Arps single segment decline or the 2-segment Modified Arps decline.
1.3 Autofit ranges
- Default: starting point for autofit - \(q_i\) is the initial production ratio, b is the rate exponent or the b-factor, and \(a_i\) is the nominal decline rate.
- Min: lower bounds of best-fit solution for \(q_i\), b and \(a_i\).
- Max: upper bound of best-fit solution for \(q_i\), b and \(a_i\).
What is the initial production ratio?
This is the ratio between the \(q_i\) that the DCA will use and the maximum production rate found in the production dataset for the well. As an example, if you have a given well where the maximum oil production in history is 1000 STB/d, having an initial production ratio between 0.8 and 1.2 means we are allowing the initial rate in DCA to be between 800 STB/d and 1200 STB/d. It is basically a generic way to specify the range for \(q_i\), for a large number of wells, and all the different phases (oil, gas and water) in one go.
1.4 Forecast
- Total Production Period: sum of historical time and forecasted time.
- Oil Cut-off rate: the rate at which to stop the forecast.
- Limiting Effective Decline Rate: specifies when the hyperbolic decline curve transitions into an exponential decline curve at a specified limiting effective decline rate.
2. Autofit
This is where you get an overview of your autofit. In this view you can filter and manually update the individual fits.
2.1 Autofit all wells
Click "Autofit all wells" to the upper left to best-fit all the wells in a given project.
2.2 What does Error (%) represent?
whitson+ facilitates regression on parameters to minimize the sum of squares of weighted residuals in the context of observed data and corresponding predictions. The objective function in is defined as:
Here, represents the RMS residual error between observed measurement and its corresponding prediction, while is the user-assigned weighting factor for that residual. Ideally, these weighting factors should be inversely proportional to the standard deviations of the residuals. Minimizing in provides the maximum likelihood estimation of the model parameters, assuming independent and normally distributed residuals.
The default weighting factors are 1. Reported is a root-mean-square (RMS) residual error () defined as:
This metric is related to but is more easily interpreted.
The residuals are calculated as relative percentages using the formula:
where is the historical production value, is the DCA value, and is the reference value for observation .
2.3 Sort, Filter and Change Streams.
- Sort: Click the header parameter to sort the well list by that parameter.
- Search: Use the search bar to search for a specific well.
- Change stream: click OIL | GAS | WATER tabs to change between the streams.
2.4 Manually Edit
You can manually check and adjust each fit by clicking the row in the table.
- Click the row. A window showing the single-well DCA fit will appear.
- Navigate to the next well with the right arrow (→) and the previous well with the left arrow (←) on your keyboard.
What does the DCA b-factor represent?
DCA fits can be tied to physical parameters as long as the fit is done while the flowing bottomhole pressure is constant (or about constant). The b-factor represent recovery efficiency. A low b-factor indicates a low recovery efficiency. A high b-factor indicates a high recovery efficiency.
2.4.1 Adjusting the autoforecast fit on a well
- Values can be changed by typing.
- The fit can also be graphically edited.
- See .gif above (it's just for example purposes).
2.5 Save Fit
Click SAVE up to the right to save all the autofits (for all the streams). These forecasts can be appended to the historical data in the type well feature.
2.6 Export
Click EXPORT up to the right to export to third-party software (for all streams).
Bulk Export DCA to ARIES
You can also use the option above to export all well declines in bulk to ARIES in one go. Typically, you want to: Export --> ARIES --> "End of Prod"
3. Summary
- Click the third step in the stepper (3. Summary).
- Click CREATE NEW CHARTS
- To quality check that the overall fits are ok, a simple "roll-up" exercise can be performed on all the data.
- That is done by adding up all the rates from all the wells producing at a given time step and plotting that over time in both rate and cumulative space.
- Here you like to visually ensure that there are no strange discontinuities when the data goes from the historical data to the forecasted part of the data.
4. Go back to the Auto-Forecast Overview?
If you want to go back to the auto-forecast overview, either 1. click the Auto-Forecast tab in the navigation panel. 2. or click the Auto-Forecast menu in the top bar.
5. Want edit Auto-Forecast well selection?
If you want to edit the well selection, either
- click the Edit Auto-Forecast tab in the Auto-Forecast Overview.
- or click the Edit Auto-Forecast menu in the top bar.
Decline Curve Analysis
Decline curve analysis (DCA) for all three flowing phases can be performed within whitson+, both automatically and manually. User-friendly, "drag-and-drop" curve-fitting functionality is also available.
1. Input
The inputs for this analysis are:
- DCA parameters (\(b\), \(a_i\), \(q_i\))
- Start (day or month)
- Total production period (years)
- Limited Decline Rate (\(a_{lim}\) or \(d_{lim}\)). More info here.
The user can decide whether to lock in the start day and the DCA parameters, or if these should change freely during the curve fitting.
2. Theory
2.1 Arps Fundamentals
The classical equation that provides that basis for this module is
where
\(q_i\) = initial flow rate
\(a_i\) = nominal decline rate at time zero
\(b\) = rate exponent
\(q(t)\) = flow rate at time t
Equation 1 is a three-parameter performance-matching equation. When the \(b\) is equal to zero, the flow rate decline is exponential, when \(b\) is 1, the rate decline is called harmonic, and when b is between zero and one, the rate decline is hyperbolic.
2.1.1 Exponential Decline (b=0)
where \(N(\infty)\) is the ultimate cumulative production.
2.1.2 Hyperbolic Decline (0 < b < 1)
2.1.3 Harmonic Decline (b = 1)
2.1.4 Beyond Hyperbolic (b > 1)
2.2 Other Segment Types
2.2.1 Linear
2.2.2 Constant
2.2.3 Semi-log
2.2.4 Power Law
2.3 DCA Summary
Type | Decline Rate | Producing Rate, q | Elapsed Time, t | Cum. Production, \(Q_{t}\) |
---|---|---|---|---|
Exponential | \(a_t = ln(\frac{q_i}{q_t})/t\) | \(q_ie^{-a_it}\) | \(ln(\frac{q_i}{q_t})/a_i\) | \(\frac{q_i-q_t}{a_i}\) |
Hyperbolic | \(\frac{a_i}{a_t} = (\frac{q_i}{q_t})^b\) | \(q_i(1+ba_it)^{-1/b}\) | \(\frac{(q_i/q_t)^b-1}{ba_i}\) | \(\frac{q_i}{a_i(1-b)}(1-\frac{q_t}{q_i}^{1-b})\) |
Harmonic | \(\frac{a_i}{a_t} = \frac{q_i}{q_t}\) | \(q_i(1+a_it)^{-1}\) | \(\frac{(q_i-q_t)}{a_iq_t}\) | \(\frac{q_i}{a_i}ln(\frac{q_i}{q_t})\) |
Linear | \(a_t = \frac{-m}{q_t}\) | \(q_i + m t\) | \(\frac{q_t - q_i}{m}\) | \(\frac{1}{m} \left[q_i(q_t - q_i) + \frac{(q_t - q_i)^2}{2}\right]\) |
Constant | \(a_t = 0\) | \(q_i\) | \(t\) | \(q_i t\) |
Semi-Log | \(a_t = ln(\frac{q_i}{q_t})/t\) | \(q_ie^{-a_it}\) | \(ln(\frac{q_i}{q_t})/a_i\) | \(\frac{q_i-q_t}{a_i}\) |
Power-Law | \(a_t = n \left(\frac{q_t}{q_i}\right)\) | \(q_i \left(\frac{q_i}{q_t}\right)^n\) | \(t = \left(\frac{q_i}{q_t}\right)^{1/n}\) | \(\frac{q_i}{1-n} \left[\left(\frac{q_i}{q_t}\right)^{\frac{1-n}{n}} - t_0^{1-n}\right], \text{ if } n \neq 1\) or \(\frac{q_i}{n} \ln\left(\frac{q_i}{q_t}\right) - q_i \ln(t_0), n = 1\) |
Type | Cumulative Production | Cumulative Production (with time substituted) |
---|---|---|
Exponential | ||
Hyperbolic | ||
Harmonic | ||
Linear | ||
Constant | ||
Semi-Log | ||
Power-Law | or | or |
2.3 Effective Decline
Nominal annual decline rate is denoted "a". The effective annual decline rate is denoted "d", with units %/yr. There are two common ways the industry converts annual decline into effective decline,
"Secant Effective" aka "Effective Hyperbolic Decline"
"Tangent Effective" aka "Effective Exponential Decline"
Both of these properties can be converted back by the following equations
Note
These numbers are mere conversions to and from \(a_{i}\) and should yield identical results irrespective of what is chosen. For most folks, it's mainly the matter of preference when it comes to what is communicated in terms of the final results.
The symbols for the different terms can vary somewhat from software to software, so here is a little overview for some handpicked softwares.
Name | whitson+ | ARIES | Harmony | PhdWin / SPEE | Val Nav |
---|---|---|---|---|---|
Initial nominal annual decline rate | \(a_{i}\) | \(A_{i}\) | \(a_{i}\) | \(D_{i}\) | Di (nom) |
Secant effective annual decline rate | \(d_{sec}\) | \(D_{b}\) | \(d_{i,sec}\) or \(d_{lim,hyp}\) | \(D_{esi}\) | Di (eff. sec) |
Tangent effective annual decline rate | \(d_{tan}\) | \(D_{h}\) | \(d_{i,tan}\) or \(d_{lim,exp}\) | \(D_{ei}\) | Di (eff. tan) |
Rate exponent | b | n | b | b | Ni |
More information about this is provided by SPEE.
2.4 Limited / Modified Decline
The limited decline rate begins when the hyperbolic decline curve transitions into an exponential decline curve at a specified limiting effective decline rate, \(d_{lim}\). The limiting effective decline rate is converted to a limiting nominal decline rate, \(a_{lim}\), and the following equations are applied:
When \(d_{tan}\) > \(d_{lim}\):
when \(d_{tan}\) <= \(d_{lim}\):
where
2.5 Plot Settings and Straight Lines
Linear, Semi-Log, and Log-Log plots provide distinct visualizations to identify decline behavior and interpret reservoir characteristics. Each plot highlights specific relationships between production rate and time or cumulative production.
Linear: Straight line on a Cartesian plot. Useful for identifying constant rate declines or trends in linear relationships. Semi-Log: Straight line on a semi-log plot. Exponential declines, often typical in fluid-flow-dominated regimes, become linear in this format. Log-Log: Straight line on a log-log plot. Power-law declines, as observed in fractured reservoirs, are best visualized in this plot.
2.6 Fetkovich
Mike Fetkovich advanced Decline Curve Analysis (DCA) by linking it directly to reservoir and physical parameters, with the b-factor playing a central role in this approach. For the b-factor to have physical significance, it must be derived while the well is producing under near-constant bottomhole pressure. The b-factor serves as an indicator of recovery efficiency, or reservoir "energy." A low b-factor signifies a low-energy system with lower recovery efficiency, while a high b-factor represents a higher-energy system with greater recovery potential. The table below outlines general guidelines for interpreting the b-factor and its relationship to reservoir characteristics.
\(b\) value | Context |
---|---|
0 | Single-phase oil, or single-phase gas at high pressure |
0.1-0.4 | Solution gas drive |
0.4-0.5 | Single-phase gas flow |
0.5-1 | Layered reservoirs |
>1 | Infinite acting flow |
2.7 Ratio Forecasting
2.7.1 Changing Forecast Type
You can change from considering the primary
2.7.2 Derivation: Why should Ratios be Power-Law?
2.7.3. How to convert Power-Law to Arps?
For further information, please refer to [1].
3. Decline Curve Analysis in whitson+
3.1 Autofit
The autofit algorithm finds the last segment of the data that is consecutively (continously) non-zero, and fits the DCA parameters to this part of the data only. A single or set of zero rates in the data is typically associated with a shut-in, in which re-initialization of the DCA fit should be performed.
3.1.1 Residual Function
The fitting routine uses the following residual function in the least square routine.
where \(r_i\) is the residual at time \(i\), \(w_i\) is the weight factor of time \(i\), \(t_0\) is the time index and q is the actual rate.
3.1.2 Default Bounds
Parameter | Lower Bound | Higher Bound | Unit |
---|---|---|---|
\(q_i\) | 40% of max observed q | 20% higher than max observed q | volume unit / time |
\(a_i\) | 0 | 10 | /yr |
\(b\) | 0 | 2 | - |
3.1.3 Adjusting default autofit bounds on DCA parameters
Fig. 1: Autofitting DCA Segments and Autofit settings
You can also adjust the settings for autofitting decline to your data by adjusting the decline parameter ranges and specifying the autofit type as the default Arps decline or the modified Arps decline.
The autofit function will automatically add a second decline segment (if it does not exist) and fit 2 segments to the data if the autofit type is modified Arps decline.
You can choose to autofit some or all of your segments at once and control the autofit on each segment by using the locks on decline parameters.
3.1.4 Weighting factors
By default, all days have a weighting factor of 1, except from the days when q = 0, in which the weighting factor is 0.
3.2 DCA on physics-based forecasts
Fig. 2: Fitting DCA to simulated history match from Numerical Model
If a physics based forecast has been established for this well via the numerical model, an option to include this data will appear as a toggle, shown above.
You can build a physics based forecast in the Numerical Model feature, where the reservoir simulation can history match the well and generate a forecast based on the numerical model set up. Once the model is simulated with a forecast included, you can bring this history match into the DCA feature to be able to build declines for the simulated history match data with the simulated forecast included, thus, allowing for physics-based, DCA forecasts.
3.3 General tips
3.3.1 Adding multiple decline segments
Fig. 3: Adding DCA Segments
You can add multiple Arps decline segments to the data, as proposed by the Modified Arps decline method (2 segments). Each of these segments can be linear incline, constant or the Arps decline type segments - controlled individually by directly entering the inputs, or adjusting the curve on the plot itself to auto-compute the decline parameters for each of them.
Adding multiple segments may be useful to model particular flow regimes, constant rate well control or refrac analysis.
3.3.3 Where does the DCA forecast append to historical data?
Fig. 4: Appending forecasts to data and EUR calculation
The end date of the forecast on the DCA plot dictates where the forecast is appended to historical data. The forecast is appended to the last producing day of the well history, or the end date of the last segment itself, whichever is earlier, and the EUR is calculated accordingly.
End date beyond the last producing day will include all the production history, followed by the forecast to calculate EUR and any end date before the last producing day will append the forecast from the set end date and follow the forecast from this point, discarding the production data after this end date.
EUR shown in the top-left box should be exactly equal to the EUR (black line) plotted on the cumulative plot (bottom) at the end of the total production period.
3.3.4 Rate-time decline vs Rate-cum decline
Fig. 5: DCA segments on rate vs cumulative plots
You can observe the decline fit in rate vs time space or rate vs cumulative space on the same plot (top).
3.3.5 Linking decline segments
You can use the 'Link to Previous Segment' checkbox for each segment to link the previous segment to this one in time for continuity. We are working on linking other elements from the previous segment (like the rates and declines at the end of the previous segment) to add to this functionality.
3.3.6 Saving the DCA fit
Fig. 6: Saving the DCA fit
You can save the current DCA fit as a saved case by using the SAVE CASE button. You can then adjust the fit as you'd like because you can always use the saved case as a checkpoint.
You can overlay the saved cases in the plot by checking the box under "Show" column for each saved case.
You can also reload the DCA parameters (and hence the fit) at any time by using the 'Load Case' option for available for every saved case by clicking the three vertical dots.
Saving a case with a name that already exists in the list of saved cases will prompt you to confirm that you want to rewrite the saved case with the new data to prevent duplicate cases.
Importing saved Type well DCA case to DCA and saved cases from DCA into Type well
If you fit DCA to the type well, made from a set of wells in the Type Well feature (DCA tab), and save that, you can import this saved DCA fit for the type well into list of saved cases in this feature. This lets you plot your type well decline over the actual well decline to compare.
Alternatively, if you saved the DCA fits for all the wells with the same name, you can import this into the Type well feature such that each of the individual well declines set up with that name is appended to the respective wells and the Type well is built for historical data + individual well declines. This removes the need for forecasting the constructed type well with DCA.
3.4 Hotkeys
The following hotkeys are available to streamline your workflow and enhance navigation within the software:
3.4.1 Decline Curve Analysis (DCA)
- Perform Autofit: Press
Alt + Enter
to automatically fit the decline curve to the available data. - Save Case: Use
Ctrl + S
to save the current case, ensuring no changes are lost.
3.4.2 Plot Navigation
- Shift Entire DCA Segment: Hold
Alt
and drag the mouse to shift the entire decline curve analysis segment on the graph. - Zoom: Click and drag the mouse over a region to zoom into specific areas of the plot for detailed analysis.
- Reset Zoom: Double-click anywhere on the plot to reset the zoom level to the default view.
3.4.3 General Navigation
- Switch Well: Use
PgUp/PgDn
to move between wells in the dataset. - Cycle to Next Phase: Press
Shift + S
to cycle between different production phases. - Change Phase: Use
Shift + O/G/W
to switch between oil, gas, and water phases.
4. Export
The results from the DCA feature can to be exported to:
- Aries
- ComboCurve
- Excel
- PHDWin
- Mosaic
- Val Nav
Note that the forecast in the export starts from the end of the uploaded historical data. Hence the best fit qi and decline parameters has to be recalculated for the end of history (as they change with time), as outlined here.
PHDWin: Uses Secant Effective or Tangent Effective Decline?
PHDWin uses Initial Tangent Effective Decline Rate.
4.1 ARIES
4.1.1 Aries Examples
The following outlines how the ARIES export works with a specific example.
For an example well with these parameters.
- qi = 51.7334 STB/d
- d_lim = 5.0000 %/yr
- b = 1.2
- d_sec = 43.7863 %/yr
- Cutoff rate = 3.0 STB/d
- forecast_time = 582.3934 months
There will be 4 cases, these are all shown assuming the parameters are for the oil phase, but the logic is equally applicable for gas and water.
Case 1 - well reaches cutoff rate before d_lim, the ARIES export would look like this:
We have a single line which shows qi, cutoff rate, b and d
DEMOWELL
PRODUCTION
CUMS 80.1645 0.4055 780.3532
START 03/2022
OIL 51.7334 3.0 B/D X MOS B/1.2 43.7863
Case 2 - well reaches cutoff rate after d_lim but before end of forecast period, the ARIES export would look like this:
We have two lines as before, only difference being that instead of putting the period of 582 months (or any number for the case), we put X
DEMOWELL
PRODUCTION
CUMS 80.1645 0.4055 780.3532
START 03/2022
OIL 51.7334 X B/D 5.0000 EXP B/1.2 43.7863
" X 3.0 B/D X MOS EXP 5.0000
Case 3 - well reaches d_lim and does not reach cut off rate until the end of the forecast period, the ARIES export would look like this:
DEMOWELL
PRODUCTION
CUMS 80.1645 0.4055 780.3532
START 03/2022
OIL 51.7334 X B/D 5.0000 EXP B/1.2 43.7863
" X X B/D 582.3934 MOS EXP 5.0000
Case 4 - d_sec is lower than d_lim:
In addition to these 3 main cases, there is an exceptional case in which d_sec < d_lim, for example if in the case presented above d_sec was 3.2561 %/yr instead 43.7863 %/yr. The export to Aries would look as follows.
Case 4a - If cutoff rate is not reached before end of forecast period:
DEMOWELL
PRODUCTION
CUMS 80.1645 0.4055 780.3532
START 03/2022
OIL 51.7334 X B/D 582.3934 MOS EXP 3.2561
Case 4b - If cutoff rate is reached before end of forecast period:
DEMOWELL
PRODUCTION
CUMS 80.1645 0.4055 780.3532
START 03/2022
OIL 51.7334 3.0 B/D X MOS EXP 3.2561
4.1.2 Aries: Secant Effective or Tangent Effective Decline?
In ARIES, the most commonly used format is 'secant effective,' also known as the 'B format.' This is what whitson+ will export to ARIES by default. However, there is also an option to use 'tangent effective,' referred to as the 'H format' in ARIES. Refer to the picture below to see the difference between these formats in ARIES.